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Abstract 

Inspired by the recent LHC Higgs data and null search results of super symmetry (SUSY), we 
scan the parameter space of the Minimal Supersymmetric Standard Model (MSSM) with relatively 
heavy sparticles (1-3 TeV). Then in the parameter space allowed by current collider experiments 



and dark matter detections, we calculate the complete one-loop SUSY QCD corrections to pp — > bbh 



at the LHC with y/s = 14 TeV and obtain the following observations: (i) For the large values of 
tan/3 and low values of tha, the SUSY QCD effects can be quite large, which, however, have been 
excluded by the latest results of LHC search for H/A — > t + t~; (ii) For modest values of tan/3 and 
mA which so far survived all experimental constraints, the SUSY QCD corrections can maximally 
reach about —9%. 

PACS numbers: 
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I. INTRODUCTION 



Very recently the ATLAS and CMS collaborations have independently reported the ob 



servation of a Higgs-like resonance with a mass about 125 GeV [1J. At the same time, the 
CDF and DO collaborations have also updated their combined results for the Higgs searches 
in bb channel, which support the LHC observation 2j. Since in the Minimal Supersymmetric 
Standard Model (MSSM) a SM-like Higgs boson is predicted with a mass below 130 GeV, 
the observation of such a 125 GeV Higgs boson supports SUSY, albeit quite restrictive on 
the parameter space of SUSY {3}. 

Meanwhile, the direct searches for SUSY particles (sparticles) have been performed at the 
LHC. Based on about 5 fb _1 luminosity, the ATLAS and CMS collaborations have reported 
null results and obtained some bounds on the sparticle masses, which is about 1 TeV for the 
gluino and first generation of squarks [4], 330 GeV for the electroweak gauginos, 180 GeV 
for the sleptons [5], 465 GeV for the stops and 480 GeV for the sbottoms [6]. These bounds 
indicate that SUSY may be heavier than expected and the sparticles may be significantly 
heavier than the electroweak scale 0,0]- 

In case that the sparticles are heavy and beyond the LHC scope of direct production, 
search for the indirect SUSY effects via loop corrections will be of great importance. Since 
the loop effects of heavy sparticles are usually small, we should look for some processes in 
which the heavy sparticles have non-decoupling loop effects. One type of such processes 
are Higgs productions at the LHC, such as the production of tH~ and hbb, in which the 
heavy sparticles have sizable non-decoupling loop effects for a small value of and a large 
value of tan/3 0, [uj] (when getting large, such effects will vanish). The reason for these 
non-decoupling effects is that the couplings in the loops are proportional to some SUSY 
mass parameters and can be enhanced by the large values of tan/3. 

In this note we focus on the production of hbb at the LHC and calculate the complete one- 
loop SUSY QCD corrections to this process. As an important Higgs production channel for 
the MSSM, this production has been studied in the literature jl0|. where the non-decoupling 
SUSY QCD effects are found to be large (reach -40% for tan/3 = 30). We revisit this 
production for the following reasons: Firstly, in the literature the SUSY QCD corrections 
to this process are calculated only partially (only the corrections to the hbb vertex have 
been considered). The complete one-loop corrections involve pentagon Feynman diagrams, 
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whose calculations are rather complicated and have not been performed. Secondly, the 
CMS collaboration has recently measured this channel and given constraints on the plane 
of tan/3 versus Since the non-decoupling SUSY QCD effects in this production is 

sensitive to the values of tan/3 and rriA, we should update the calculations by considering 
such new constraints. Moreover, other experimental constraints, such as the dark matter 
direct detection limits and the SM-like Higgs boson mass around 125 GeV, are also rather 
restrictive and should be considered. In this work, we consdier all current experimental 
constraints to scan the MSSM parameter space and then in the allowed parameter space we 
calculate the process pp — > bbh with the complete one-loop SUSY QCD corrections. 

The paper is organized as follows. In Sec. II. we describe the calculations for the process 
pp — » bbh. In Sec. Ill we show numerical results. Finally, we draw the conclusions in Sec. 
IV. 



II. THE DESCRIPTION OF CALCULATIONS 

In the MSSM the lighter CP-even Higgs mass (m^) is smaller than Mz at tree level but 
receives large corrections at the loop level. The leading part of the corrections is from the 
stop sector and can be expressed as [kJ 

where X t = A t —fi cot f3 is the stop mixing parameter. We see that a large stop mass or a large 
stop mixing parameter is needed to increase to 125 GeV. In our calculations we consider 
the collider constraints on the MSSM Higgs sector, using the packages FeynHiggs2.8.6 [l^] 



and HiggsBounds-3.8.0 15] to calculate the observables in the Higgs sector and require them 
to satisfy the constraints from the LEP, Tevatron and LHC. 

The SUSY QCD corrections to hbb production at the LHC involve the sbottoms and 



gluino in the loops. The sbottom mass matrix takes the form 



16] 
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A b — fi tan {3 
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(4) 
(5) 



where mj~ and m 2 ^ are respectively the soft-breaking mass parameters for the left-handed 
squark doublet Q and the right-handed down squark D, A b is the sbottom soft-breaking 
trilinear coupling and \x is the SUSY-preserving bilinear coupling of the two Higgs doublets 
in the superpotential. This mass matrix can be diagonalized by a unitary transformation 
which rotates the weak eigenstates Slr to the mass eigenstates b\ 2, 
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with the sbottom masses m~ bi 2 and the mixing angle 6~ b determined by 
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We produce the one-loop amplitudes with FeynArts-3.5 \v\ and use the FormCalc-6.1 [l^ 
to simplify them and express the loop functions defined in [191 ] . The numerical calculations 
are performed by using LoopTools-2.2 [20|]. In order to preserve supersymmetry, we adopt 
the constrained differential renormalization (CDR) 21[ to regulate the ultraviolet divergence 
(UV) in the self-energy and vertex corrections, which is equivalent to the dimensional reduc- 



tion method at one- loop level 22]. In Fig.l we display the representative pentagon Feynman 
diagrams for the SUSY QCD corrections in the subprocesses gg — » bbh. Due to no massless 
particles in the loop, all the Feynman diagrams with the gluino and sbottoms in the loops 
are infrared (IR) finite. 



In our calculations, we assume a common SUSY mass Msusy defined by M t 



SUSY 



M, 



Q 



Mfj = = Mg — At — Af, — jj,. We fix slepton mass parameters Mf = = A T = 3 TeV 
and scan the following MSSM parameter regions: 

5<tan/3< 60, 90 GeV < M A < 350 GeV, 1 TeV < M SU sy < 3 TeV (9) 
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FIG. 1: The pentagon diagrams for SUSY QCD corrections to the subprocess of gg — > bbh at the 
LHC. 

In our scan we consider the following constraints on the parameter space: (i) We require that 
the SM-like Higgs is in the region of 120 GeV< < 130 GeV; (ii) For the constraints from 
flavor physics and electroweak precision tata, we checked hy using the package susy.flavor 
v2.0 [23] that they are safely satisfied because we assume relatively heavy sparticles. (iii) We 
consider the dark matter constraints from the WMAP relic density and the direct detection 
results by using the package MicrOmega v2.4 [24]. 



III. NUMERICAL RESULTS 



Since the b-quark Yukawa coupling may receive large radiative corrections in the MSSM, 
we use the running b-quark mass (m^ R ) and use the method induced in 26] to absorb the 
MSSM corrections into the effective b-quark Yukawa couplings. But for the b-quark in the 
final state, we take the pole mass to assure the correct on-shell behavior. 



In our numerical calculations, we take the input parameters of the SM as 



27] 



m t = 172 GeV, mf s (mf s ) = 4.19 GeV, m z = 91.1876 GeV, a(m z ) = 1/127.918 



Here m^ IS (m MS ) is the QCD-MS 1 bottom-quark mass, which is related to m^ R as 

" 2 -(73-3 % )] 
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where rif is the number of the active quark flavors. For the strong coupling constant 
we take its 2-loop evolution with QCD parameter A nf=5 = 226 MeV and get a s (mz) = 0.118. 
We use CTEQ6L1 and CTEQ6M [28] parton distribution functions (PDF) for the SM tree 
level and SUSY QCD one-loop level computations, respectively. The renormalization scale 
[ir and factorization scale \if are chosen to be /jlr — Hf — rnz- We numerically checked 
that all the UV divergence in the loop corrections canceled. 
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FIG. 2: The scatter plots of the samples satisfying constraints (i)-(iii), projected on the planes of 
tan /3 and versus rriA- The curves denote the bounds from the measurements of pp — > H/A — > 
t + t~ and pp — > bbh(h — > bb) at the LHC (the region above each curve is excluded), The triangles 
(red) represent the samples survied all the constraints. 

In Fig.2 we project the survived samples satisfying the constraints (i)-(iii) on the planes 
of tan/3 and versus m^. It can be seen that they spread in three different regions: 
the non-decoupling regime (m^ < 200GeV, tan/3 > 20), the decoupling regime with a 
large tan/3 {m A > 200GeV, tan/3 > 40) and the decoupling regime with a small tan/3 
(rriA < 200GeV, tan/3 < 15). In this figure we also show the bounds from the Higgs searches 
pp — > H/A — > t + t~ [l2j and pp — > bbh at the LHC [ll|. We can see that the constraint 



from the t~t + channel is much stronger than pp — > bbh and two regimes with a large tan /3 
have been excluded. Only the decoupling regime with a small tan /3 is still allowed, where 
the SM-like Higgs can get a correct mh ~ 125 ± 2 due to the contribution of heavy stops. 



-20 



^ -40 

b 

to 



-60 



-80 



-i 1 1 1 1 1 1 1 1 r 

- r »4att»» 



-1 1 1 1 1 1 1 1 1 1 1 r 

-- .v- • • • » % -r9 



•• • 



4** 



M A (GeV) 



tan/3 



i — i — 1 — i — < — i — < — i — < — i — i — < — i — 1 — i — < — i — < — i — < — i — <- 

100 150 200 250 300 350 10 20 30 40 50 60 



FIG. 3: Same as Fig. 2, but showing the one-loop SUSY QCD corrections. 



In Fig. 3 we present the one-loop SUSY QCD corrections to the process pp — > bbh for the 
samples displayed in Fig. 2. We can see that the SUSY QCD corrections will be significant 
for the samples which have a large tan/3 and a low value of uia- This can be understood 
by the contribution to the effective b-quark Yukawa coupling after integration of the heavy 



sparticles, which is Sy 



hbb 



ja s mj"sina/ Mgfi 
6nmw cos /3 \ M 



SUSY 



) [tan /3+cot a] [8(. Since we assume Msusy 



/i, the Yukawa coupling will be independent of the sparticle masses and be greatly 



enhanced by a large tan (3. However, it should be noted that these samples will lead to the 
excess of the production rate of pp — > H/A — > t + t~ and thus have been excluded by the 
current measurements. With the increase of and the decrease of tan {3, the corrections 
drop rapidly and approach zero in the decoupling limit. The main reason is that 5y hb i can 
be heavily reduced by the cancellation between tan/3 and cos a, which have a relation as 
cot a ~ — tan/3 — 2m| tan /3 cos 2 f3/m\ for a large m^. For the samples which survived all 
the constraints, the SUSY QCD corrections can only reach about —9% at the LHC with 
^/i = 14 TeV. Detecting such a size of SUSY QCD effects may be challenging in the future 
measurement of the process pp — > bbh [29| . 
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IV. CONCLUSION 



In this work, we calculated the complete one-loop SUSY QCD corrections to the process 
pp — > bbh at the LHC with a/s = 14 TeV. We found that the large SUSY QCD corrections in 
the non-decoupling regime with a large tan /3 and a low uia has been excluded by the latest 
results of LHC non-standard Higgs searches. For the survived decoupling regime which have 
modest values of tan/3 and uia, the SUSY QCD corrections can maximally reach —9%. 
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